Nonleptonic Q decays and the Skyrme model 
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Nonleptonic Q~ decay branching ratios are estimated by 
means of the QCD enhanced effective weak Hamiltonian sup- 
plemented by the SU(3) Skyrme model used to estimate the 
nonperturbative matrix elements. The model has only one 
free parameter, namely the Skyrme charge e, which is fixed 
through the experimental values of the octet-decuplet mass 
splitting A and the axial coupling constant </a- The whole 
scheme is equivalent to that which works well for the non- 
leptonic hyperon decays. The ratios of calculated amplitudes 
are in agreement with experiment. However, the absolute val- 
ues are about twice too large if short-distance corrections and 
only ground intermediate states are included. 

PACS number(s): 12.39. Dc, 12.39. Fe, 13.30.Eg 



Recently both s- and p-wave nonleptonic hyperon de- 
cay amplitudes were quite successfully reproduced by the 
SU(3) extended Skyrme model with the QCD enhanced 
effective weak Hamiltonian [Q. The decay amplitudes 
were described through the current-algebra commuta- 
tor, the ground-state baryon pole terms and factorizable 
contributions. The nonperturbative quantities, i.e. the 
baryon 4-quark operator matrix elements, were estimated 
using the SU(3) Skyrme model. For the s-wave hyperon 
decay amplitudes correct relative signs and absolute mag- 
nitudes were obtained. For the p-waves all relative signs 
were correct, with their relative magnitudes roughly fol- 
lowing the experimental data. As far as the absolute 
magnitudes are concerned, the poorest agreement be- 
tween theory and experiment was within a factor of two. 
One is thus faced with an obvious question: could an 
analogous approach work equally well for the f2 _ nonlep- 
tonic weak decays? Such a question should be considered 
in connection with the accurate measurements of the Q~ 
decay branching ratios. The experimental value for the 
Q~ mean life time is ||: 

T CX pcrimcnt(^~) = (82.1 ± 1.1) pS 

from which the values of the p-wave amplitudes given in 
Table | have been extracted. 

In the Skyrme model, baryons emerge as soliton config- 
urations of the field U of pseudoscalar mesons ||[| . Ex- 
tension of the model to the strange sector |5|-|8| is done by 
an isospin embedding of the static hedgehog Ansatz into 



an SU(3) matrix which is a subject of a time dependent 
rotation 



U(r,t)=A(t) 



exp(if-nF(r)) 







by a collective coordinate matrix A(t) £ SU(3) 
which defines the generalized velocities A'(t)A(t) — 
| Yl a =i A Q a Q and the profile function F(r) is interpreted 
as a chiral angle that parametrizes the soliton. The col- 
lective coordinates a a are canonically quantized to gen- 
erate the states that possess the quantum numbers of 
the physical strange baryons. In order to account for a 
non-zero strange quark mass the appropriate chiral sym- 
metry breaking (SB) terms should be included. In this 
work, however, following Q we shall neglect the SB ef- 
fects since they are much smaller than the uncertainties 
coming from other sources (like for example the values of 
the Ci coefficients in the effective Hamiltonian, or factor- 
ization. Moreover the SU(3) symmetry is used in course 
of the calculations of the decay amplitudes, so for con- 
sistency the nonperturbative matrix elements have to be 
evaluated in the same approximation scheme. 

Our goal is to see whether the effective weak Hamilto- 
nian and the SU(3) Skyrme model are able to predict the 
nonleptonic f2 _ decay amplitudes (CIk) ■ ^ — ► Aif~, 



(fT) : Q- 



and (f2 ) : £1 



77-^0 



following 



the method of Refs. MM. To this end we shall employ 
the Standard Model effective Hamiltonian and the mini- 
mal number of couplings concept of the Skyrme model to 
estimate the nonperturbative matrix elements of the 4- 
quark operators jjlfioj and the axial current form- factor 
for decuplet-octet transition. This approach uses only 
one free parameter, i.e. the Skyrme charge e. In order 
to avoid the unnecessary numerical burden, throughout 
this report we use the arctan Ansatz for the Skyrme pro- 
file function F(r) J7pl|], which allows to calculate the 
pertinent overlap integrals analytically with accuracy of 
the order of a few % with respect to the exact numerical 
results. 

It is well known that the nonleptonic weak decays of 
baryons can be reasonably well described in the frame- 
work of the Standard Model pjl2[|. The starting point 
in such an analysis is the effective weak Hamiltonian in 
the form of the current ® current interaction, enhanced 
by Quantum Chromodynamics (QCD), i.e. obtained by 



1 



integrating out the heavy-quark and the IT-boson fields, 



H°"(AS = 1) = V2G F V* d V u , 



4 
i=l 



v: d v u . 



(2) 



are the 



where Gf is the Fermi constant and 
Cabbibo-Kobayashi-Maskawa matrix elements. The co- 
efficients Ci are the well known Wilson coefficients [ jl2| , 
most recently evaluated in Ref . |l3| and the Oi arc the 
familiar 4-quark operators fl2] , |l3|| . For the purpose of 
this work we neglect the so called Penguin operators 



since their contributions are proven to be small |12|,i3| . 
We are using the Wilson coefficients from Ref. [12|: 
ci = -1.90-0.61C ,c 2 = 0.14+0.020C ,c 3 = c 4 /5 ,c 4 = 
0.49 + 0.005C, with C = V t * d V ts /V* d V us . Without QCD- 
short distance corrections, the Wilson coefficients would 
have the following values: c\ — — 1, C2 = 1/5, C3 = 
2/15, C4 = 2/3. In this paper we simply consider both 
possibilities and compare the resulting amplitudes. 

The techniques used to describe nonleptonic de- 
cays (3/2 + — > 1/2+ + 0~ reactions involve only p- and d- 
waves) are known as a modified current-algebra (CA) ap- 
proach. The general form of the decay amplitude reads: 

<n(q)B'(p')\H e Jf\B(p)> = U(p')[B + l5 C}q»W,(p). (3) 

Here U (p') denotes regular spinor while and W M (p) is the 
Rarita-Schwinger spinor. The parity-conserving ampli- 
tudes B correspond to the p-wave and parity-violating 
amplitudes C correspond to the d-wave fl~ decays, re- 
spectively. 

The decay probability T(3/2+ -> 1/2+ + 0") reads: 



,'|3 



r = 



[(E' + m f )\B\ 2 + (E'-m f )\C\ 2 ], 



,'|2 



[(» 



+ ml)/2m n ] 2 -m\ 



E' = (m Q + mf - mj)/2mn. 



(4) 



Here rrif denotes the final baryon mass and m$ is 
the mass of the emitted meson. Obviously the parity- 
violating amplitude C (d-wave) is multiplied in T by an 
unfavorable factor [E' — rrif). In our calculation frame- 
work, the amplitude C can receive contributions only 
from the pole diagrams. As these pole diagrams con- 
tain negative-parity 1/2 - and 3/2~ baryon resonances, 
the sum of baryon masses instead of the difference ap- 
pears in denominator. Parity-violating amplitudes are 
thus suppressed by a large factor. If the vertices in both 
B and C pole contributions are of the same order of mag- 
nitude, then the decay of £l~ is almost parity-conserving. 
This conclusion is the same as the one drawn in Ref. [[14j . 

We calculate the parity-conserving p-wave tt~ decay 
amplitudes B using the so-called tree-diagram approx- 
imation at the particle level. This means that in this 
work we take into account all possible tree diagrams in- 
volving baryons and mesons, i.e. factorizable and pole 
diagrams. 



All B amplitudes receive contributions from the pole 
diagrams, as shown in Fig.l, and they are as follows: 



9Q-K-B° a A3o 9S*~ K~ AO'S— Q- 



m^o — m\ m n- 

m n - - m H »- 
ff5-s*-7i- oQ E:*-n- 
m n - - m H «- 



(5) 




FIG. 1. Pole diagrams. The double lines are 3/2 + reso- 
nances. The full lines are hyperons and the dashed lines are 
mesons. Weak and strong vertices are indicated. 

The decuplet-octet-meson strong coupling constants 
determined from the experimental value of the decay rate 
A ++ — > pn + and by using SU(3) relations are as follows: 



9A++7T + P 



9~.°n-K- 



-9a- 



15.75GeW 



(6) 



Note that the strong couplings g^Nn and gNN-K are re- 
lated by gANir = § gNN-K, which follows from the 1/N C ex- 
pansion without other assumptions [0 . This relation is in 
excellent agreement with experiment and with the evalu- 
ation of gNNn and gANn in the Skyrme model Note 
also that gNN* and g^ N are related by the Goldberger- 
Treiman relation, which is also true in the Skyrme model. 
So altogether we can reproduce the experimental value 
of gANir if we use the Skyrme model parameters which 
reproduce g^ N ■ However, consistently following the ap- 
proach of Ref. , we are using the Skyrme model to esti- 
mate only the unknown matrix elements. Therefore, it is 
natural to use the experimental value gA++a+p — 15.75 



GeV -1 in the calculations of baryon poles (|5[) 

The baryon-pole B-p amplitudes contain weak matrix 
elements defined as a,BB' = G < B'\ J2i=i c iOf c \B >, 
were the constant Q — \f2GFV* d V us . The important part 
of clbb' are 4-quark operator matrix elements, which are 
nonperturbative quantities. This is the first point of this 
work at which the Skyrme model is used. 

The factorizable contributions to p-waves (only for 
(f2Z ) amplitudes), are calculated by inserting the vac- 
uum states; it is therefore a factorized product of two 
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current matrix elements, where the octet-decuplet ma- 
trix element of the axial-vector current reads: 

< E(p')\A»\n-(j>) >=gT{q 2 )U B {j>')W^). 

Summing over all factorizable contributions gives the fol- 
lowing expressions for the amplitudes: 

B s (fT) = -^jff*9A n ~ [ci - 2(c 2 + c 3 + c 4 )] , 

B s (%) = \0UgT a ' [a - 2(c 2 + c 3 - 2c 4 )] . (7) 

The g^ 1 represents the form-factor of the spatial compo- 
nent of the axial-vector current between and f2 
states. This is the second point of this paper where 
the Skyrme model enters the calculation. Note that the 
| Bs(Qq ) I <C \Bs(QZ)\ due to the helicity suppression 
enhanced by the QCD corrections. 
Total theoretical amplitudes are: 



e w 4 was successfully adjusted to the mass difference of 
the low-lying 1/2+ and 3/2+ baryons (Table 2.1 of [§). 
This value of e was next employed to evaluate the static 
properties of baryons. 

For the evaluation of the nonlcptonic f2~ decays, the 
most important baryon static properties is the octet- 
decuplet mass splitting A. Other important quantity 
is the axial coupling constant g\ . We compute these 
quantities by using the arctan Ansatz in the SU(3) ex- 
tension of the Skyrme Lagrangian ©. By fixing A and 
g v ^ to their experimental values, we obtain e = 4.21 and 
e = 3.41, respectively. In further calculations of the f2~ 
decay amplitudes, we use the mean value e = 3.81, like 
in the case of the nonleptonic hyperon decays (!]]. As in 
the latter case this introduces approximately 15 % un- 
certainty of the decay amplitudes which are dominated 
by & SK (see (jl3| ) below) which scales like 1/e. 

We proceed with the computation of the axial-current 
form-factor g^ n in the Skyrme model using the arctan 
Ansatz: 



£thcory(<7 2 ) = B v (q 2 )+B s (q 2 ), 



(8) 



were the relative signs between pole and factorizable con- 
tributions are determined via SU(3) and the generalized 
Goldberger-Treiman relations. 

In order to estimate the matrix elements entering 
Eqs.(||) and J7), we take the SU(3) extended Skyrme La- 
grangian pip pi: 



C = C a +C 



Sk 



Csb + Cwz, 



(9) 



where L a , Csk, £-sb, and C>wz denote the cr-model, 
Skyrme, symmetry breaking (SB), and Wess-Zumino 
(WZ) terms, respectively. Their explicit form can be 
found in Refs. pjl5]]. In this work we will use the SU(3) 
extended Lagrangian (^) in the limit of fx = fir- Then 
our new set of parameters i, /3', and 8', determined 
from the masses and decay constants of the pseudoscalar 
mesons g, is: x = 2m 2 K /ml - 1, $ = 0, 8' = m£/*/4. 
Owing to the presence of the 8' term in classical mass E c \ , 
the dimcnsionless size of soliton x' becomes a function 
of e, fa, and 8': 



r' 2 
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(10) 



Note that x is responsible for the baryon mass split- 
tings and the admixture of higher representations in the 
baryon wave functions. However, as explained above, in 
this work we use the SU(3) symmetric baryon wave func- 
tions in the spirit of the perturbative approach to SB. 
Indeed, we have shown in Ref. Q, that the SB effects 
through the weak operators in the decay amplitudes are 
very small. 

Our fitting procedure is as follows. Since the coupling 
/tt is equal to its experimental value, the only remain- 
ing free parameter is the Skyrme charge e. The value 



9 T n ~ 



9a 



ff ! n K) = -x/l5 4 n K), (11) 



were g A n {x' ) is given in eq.(10) of Ref . Q. 

Let us recapitulate results of Ref. E| where we have 
calculated the matrix element of the product of two (V — 
A) currents between the octet states using the Clebsch- 
Gordon decomposition |10|Jl5|| : 



(12) 



<B 2 \OW\B 1 >=<^ SK xY^C R 



where § SK is a dynamical constant and Cr denote the 
pertinent sum of the SU(3) Clebsch-Gordon coefficients 
in the intermediate representation R. The same holds 
for the WZ current. The total matrix element is then 
simply a sum < o\ SK ^ + > 5 with i — 1, ... ,4. 

The quantities <E> are given by the overlap integrals of 
the profile function in Ref. Q. 

For the 0\ operator, R — 8 QiS or 27. The matrix ele- 
ments read: 



< A° /2 |Oi|2° /2 >= (18.14$ SK - 20.41$ 



WZ' 



10" 



= (4.86|sif - 0.09|M/z)10- 3 GeV 3 , (13) 
< 3*7 2 |6i|f2- /2 >= (-10.60$ 5A ' + 12.89$ wz )10- 3 
= (-2.84| SJC + 0.06|M/z)10- 3 GeV 3 . (14) 

Note that we have omitted the static kaon fluctuations 
in all of our computations. 

Our main goal was to learn how the approach of Refs. 
jl],[)| in which the Skyrme model is used to estimate 
the unknown nonperturbative matrix elements applies 
to the 0~ nonleptonic decays. We have systematically 
presented all possible tree diagrams, namely factorizable 
contributions, octet diagram (only for (fix) amplitude) 
and decuplet pole diagrams. Numerical results presented 



3 



in Table | are very encouraging. They are in satisfac- 
tory quantitative agreement with experimental data. In 
Table | both factorizable (Bs(q 2 )) and pole-amplitudes 
(B-p(q Y) are displayed. Comparison of the total ampli- 
tudes (g) #theory(<Z 2 ) with experiment shows the follow- 
ing: 

(a) Dynamics based on the pole-diagrams supported 
by the Skyrme model leads to very good results for rel- 
ative importance of various decay modes. However, all 
amplitudes are too large by about a factor of two. 

(b) The decay amplitude £Ik , which does not contain 
factorizable contributions, has the largest pole contribu- 
tion, since both pole diagrams are of approximately the 
same size and contribute constructively. 

(c) Short-distance corrections to the effective weak 
Hamiltonian are proved to be important |^|,[l| . 

(d) In fact, too large values of the QCD-enhanced 
Bj>(q 2 ) amplitudes are actually a welcome feature. 
Namely, if one takes into account the higher S(3/2 + ) 
resonances then, the poles (||) would flip the signs, since 
their masses are larger than the f2~~ mass. The assump- 
tion that the vertices would be more or less the same 
like for the ground-state poles leads to internal cancella- 
tions between pole-diagram contributions ■ Taking 
into account (a) - (c) and Ref. Q we expect that this 
dynamical scheme would produce better agreement with 
experiment. 

(e) We have found the octet dominance i.e. the 27- 
cont animations are small. 

(f) We also find the dominance of the Skyrme La- 
grangian currents over the WZ current. For e ~ 4, the 
WZ contribution to ((l|Jl|) is below 3%. 

(g) Like in the nonleptonic hyperon decays, the factor- 
izable contributions, in the Skyrme model approach, turn 
out to be not important for the nonleptonic decays. 

(h) Finally, we have found that the pole terms and the 
factorizable contributions have opposite signs. 

We conclude that in general our approach provides the 
good description of the f2~ decays. Obviously, not all 
details are under full control (e.g. m s corrections are ne- 
glected); nevertheless, it seems the QCD-corrected weak 
Hamiltonian H^ff, together with the inclusion of other 
possible types of contribution to the total amplitudes 
(-KT-poles and/or factorization, higher baryon-poles, etc.) 
supplemented by the Skyrme model, lead to a correct 
answer. 

We hope that the present calculation, taken together 
with the analogous calculation of the nonleptonic hy- 
peron decay amplitudes will contribute to the un- 
derstanding of the nonleptonic hyperon interactions. It 
is also a test of the application of the Skyrme model 
to the evaluation of the nonperturbative quantities like 
axial coupling constants and the dimension six operator 
matrix elements between the different baryon states. 
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TABLE I. The p-wave (B) nonleptonic £l~ decay ampli- 
tudes. Choices (off, on) correspond to the amplitudes without 
and with inclusion of short-distance corrections, respectively. 
For Q.k (^o"_) amplitude q 2 = rn 2 K (rn 2 T ). 



Amplitude (lO^GeW 1 ) 


(si K ) (fii) ov) 


< 2+ V) 


off 2.67 
on 5.51 


4 3/2+ V) 


off 1.57 1.28 0.91 
on 3.14 2.56 1.81 


Bs{q 2 ) 


off -0.27 -0.06 
on -0.30 0.03 


Bthcory(g 2 ) 


off 4.24 1.01 0.85 
on 8.65 2.26 1.84 


ScxpUO^GeV- 1 ) § 


4.02 1.35 0.80 
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